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METHODE POUR DEFORMER GRADUELLEMENT 
HETEROGENE TEL QU'UNE ZONE SOUTERRAINE. 

(57) . Methods pour defomier graduellement Ies represen- 
tations ou realisations. g6n6r6es par simulation sequentlel- 
le, d'un rriodfeie stochastique non n^cessalrement gausslen 
d'une grandeur physique z dans un milieu h6t§rog6ne 
maill6, afin de Ies ajuster a un ensemble de donnees relati- 
ves k la structure ou l'6tat du milieu qui sont collectees par 
des mesures et observations pr^alables. 

- Elle comporte essentiellement rappllcation d'un algo- 
rithme de deformation graduelle d'un module stochastique 
h un vecteur gausslen h N variables mutuellement Ind6pen- 
dantes, qui est relie k un vecteur uniforms ^ N variables uni- 
fonmes mutuellement ind6pendantes par la fonotlon de 
repartition gausslenne de facon a ddfinir des realisations du 
vecteur unlforme, et {'utilisation de ces realisations pour ge- 
nerer des representations de cette grandeur physique z, 
que I'on cale par rapport aux donnees. 

- Applications par exemple pour visualiser la configura- 
tion statlstique d'une grandeur: la permeabillte d'un glse- 
ment souterrain, la pollution atmospherique, etc. 
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La pr^sente invention a pour objet une methode pour deformer 
graduellement des representations ou realisations, generes par simulation 
sequentielle, d'un moddle stochastique non necessairement gaussien d'un 
milieu heterogene, en se basant sur un algorithme de deformation graduelle de 
5 modeles stochastiques gaussiens. 

La methode selon T invention trouve des applications dans la 
modelisation de zones souterraines oii il s'agit de generer des representations 
montrant comment est distribuee une certaine grandeur physique dans une 
zone du sous-sol (la permeabilite z par exemple), qui soient compatibles au 
10 mieux avec des donn^es observdes ou mesurdes : donnees geologiques, 
enregistrements sismiques, mesures obtenues dans des puits notumment 
mesures des variations au cours du temps de la pression et du debit de fluides 
issus d' un gisement etc- 

Etat de la techmque 

15 Par la demande de brevet FR 98/09.018. on connatt une mediode pour 

deformer graduellement un modele stochastique (de type gaussien ou 
apparente) d'un milieu heterogene tel qu'une zone souterraine, contraint par 
une ensemble de param^tres relatifs a la structure du milieu. Cette methode 
comporte le tirage d'un nombre p Cp=2 par exemple) de realisations (ou 

20 representations) independantes du modele ou d'une partie au moins du modele 
choisi du milieu, parmi T ensemble de toutes les realisations possibles et une 
ou plusieurs etapes itdratives de deformation graduelle du modele en 
effectuant une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes enue elles puis des realisations composites 

25 successivement obtenues eventuellemeni avec de nouveaux tirages, etc., les 
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coefficients de cette combinaison etant tels que la somme de leurs carres est 
egaleal. 

Les modeles gaussiens ou apparentes cpnviennent pour modeliser des 
champs de grandeurs continues et, de ce fait, sont mal adaptes a la 
5 modelisation de zones traversdes par des r^seaux de fractures ou des systdmes 
de chenaux, par exemple. 

On connait par ailleurs les algoriihmes dits de simulation sequentielle, 
les plus largement utilises parmi les algorithmes de simulation geostatistique. 
Bien qu'ils conviennent particulierement pour la simulation des modeles 
10 gaussiens, il n'impliquent pas en principe une limitation h ce type de modele. 

On forme une representation geostatistique d*une zone souterraine par 
exemple en la subdivisant par un reseau Jl N mailles et en determinant un 

vecteur aleatoire a N dimensions Z = (Zi, Zn) correspondant au mieux i 

des mesures ou observations obtenues sur la zone. Comme le montre par 

15 exemple Johnson, M.E. ; in « Multivariate Statistical Simulation » ; Wiley & 
Sons, New York, 1987, cette approche reduit le probleme de la creation d'un 
vecteur a N dimensions i une sdrie de N problemes mono-dimensionnels, Un 
tel vecteur aldatoire n'est ni necessairement multi-gaussien ni stationnaire. La 
simulation sdquentielle de Z implique d'abord la definition d'un ordre selon 

20 lequel les N elements (Zi, Z^ Zn) du vecteur Z sont generes Tun apres 

r autre. Sans rien perdre de generalite, on suppose que les N elements de Z 
sont gdneres en sequence depuis Z\ jusqu'clZN. Pour tirer une valeur de chaque 
element Zi, (i = 1, X), il faut realiser les operations suivantes : 

a) construire la distribution de Ziconditionnee par (Zi, Z:.. Z.-i) 
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- P(Z,- < IZ,.Z^ Z,.j) ; et 

b) tirer une valeur de Zj de la distribution Fc(Zi). 

Dans la pratique de la geos|atistique, la simulation s^quentielle est 
utilisee frequemment pour engendrer des vecteurs multi-gaussiens et des 
5 vecteurs d'indicateur non-gaussiens. La tache principale de lai simulation 
sequentielle est la determination des^distributions conditionnelles Fc(Zi) (i = 1, 
NT). Des algorithmes et logiciels pour estimer ces distributions sont decrits 
par exemple dans : 

- Deutsch, C.V. et al, « GSLIB (Geostatistical Software Library) and User's 
10 Guide » ; Oxford University Press, New York, Oxford 1 992. 

En ce qui conceme le tirage des valeurs a partir de la distribution Fc(Zi), 
il existe aussi un vaste ensemble d' algorithmes connus. 

On considSre la mdthode de distribution inverse par laquelle on obtient 
une realisation de Zr. = F/*(w.), ou Ui est tire d'une disu-ibution uniforme 
15 entre 0 et 1. II s'ensuit qu'une realisation du vecteur Z correspond a une 
realisation du vecteur U dont les elements Ui, U2, ....,UNt sont mutuellement 
independants et uniformement distribues entre 0 et 1 . 

On s'aperQoit qu'une simulation sdquentielle est une operation S qui 
transforme un vecteur uniforme U= (Ui, U:, ...mL'n ) en un vecteur structure Z 
20 = (Zi, Z:,.... Zn) : 

Z = S(U) (1) 

Le probleme que represente la contrainte d'un vecteur Z a differents 
types de donnees peut etre resolu en contruignant les distributions 
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conditionnelles Fc(Zi) (i = 1, N) et/ou le vecteur uniforme U= (Ui, U:, 
..,„Un)- 

Des travaux recents sur ralgorithme sequentiel ont ete concentres sur 
I 'amelioration de Testimation des distributions conditionnelles F^:(Zi) par 
5 integration de donnees geologiques et de donnees sismiques. 

On peut citer par exemple le papier de Zhu, H. et al ; « Formatting and 
Integrating Soft Data : Stochastic Imaging via the Markov-Bayes Algorithm » 
in Soares, A., Ed. Geostatistics Troia 92, vol. I : Kluwer Acad. Publ., 
Dordrecht, The Netherlands, p.1-12, 1993. 

10 Cependant, cette approche ne peut ^tre etendue h, T integration de 

donnees non Uneaires comme les pressions issues des tests de puits et les 
historiques de production, sauf h imposer une severe linearisation. En outre, 
comme toute combinaison de vecteurs uniformes U ne donne pas un vecteur 
uniforme, la methode de deformation graduelle d'un modele stochastique 

15 developpde dans le brevet d6j^ cite, ne peut etre directement appliquee dans le 
cadre de la technique sdquentieile rappelee ci-dessus. 

Aussi, la methode selon T invention permet de rendre compatible les 
deux approches i.e. d'etendre le formalisme developpe dans le brevet 
prdcddemment cite, h la deformation graduelle de realisations, gdn^rees par 
20 simulation sequentielle, d'un modele stochastique non necessairement 
gaussien. 

Definition de la methode 

La methode permet de deformer graduellement une reprdsentation ou 
realisation, gendrde par simulation sequentielle, d'un modele stochastique non 
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necessairement gaussien d'une grandeur physique z dans un milieu h^terog&ne 
tel qu'une zone souterraine, afin de la contraindre h, un ensemble de donnees 
collectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa structure. 

5 Elle est caracterisee en ce qu'elle comporte T application d'un 

algorithme de deformation graduelle d'un modele stochastique un vecteur 
gaussien (Y) ayant un nombre N de variables mutuellement independantes qui 
est relie ^ un vecteur uniforme (U) h N variables uniformes mutuellement 
independantes par la fonction de repartition gaussienne (G), de fagon a defmir 

10 une chaine de realisations u(t) du vecteur (U), et T utilisation de ces 
realisations u(t) pour generer des realisations z(t) de cette grandeur physique 
que Ton cale par rapport aux donnees (non lineaires). 

Suivant un premier mode de mise en ceuvre, on defmit la cbaine de 
realisations u(t) du vecteur uniforme (V) h partir d'une combinaison lineaire 
15 de realisations du vecteur gaussien (Y) dont les coefficients de combinaison 
sont tels que la somme de leurs carr^s est ^gale k T unite. 

Suivant un autre mode de mise en osuvre, on realise une deformation 
graduelle d'un nombre n de parties du module repr^sentatif du milieu 
heterogene tout en preservant la continuite entre ces n parties du modele. par 
20 une subdivision du vecteur uniforme (U) en un nombre n de sous-vecteurs 
mutuellement independants. 

D*autres caracteristiques et avantages de la methode selon Tinvention, 
apparaitront a la lecture de la description ci-apres d'un exemple non limitatif 
de realisation, en se referant aux dessins annexes oil : 
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- la Fig.l montre la couche mediane d'une realisation d'un modele de facias 
servant de reference, generic par simulation sequentielle d* indicatrices ; 

- la Fig,2 montre la variation en fonction du temps de la pression obtenue 
dans un test de puits du modele de reference ; 

5 - les Fig.3A i 3E montrent respectivement cinq realisations initiales de la 
couche mediane d'une zone reservoir, contraintes seulement par les facids le 
long du puit ; 

- les Fig. 4A a 4E montrent respectivement pour ces cinq realisations, les 
courbes des pressions au puits dans le modele de reference compares i 

10 celles correspondant aux modeles initiaux ; 

- les Fig.SA h 5E montrent respectivement cinq realisations de la couche 
mddiane du modele de facies conditionnees aux facids le long du puits et 
calees par rapport £l la courbe de pression obtenue par test de puits ; 

les Fig. 6A a 6E montrent respectivement pour les cinq realisations, les 
15 courbes des pressions au puits dans le modele de reference compardes h 
celles correspondant aux modeles cales ; 

- les Fig.7A h 7E montrent respectivement comment les fonctions objectif 
correspondant respectivement a ces cinq exemples, varient avec le nombre 
d' iterations ; 

20 - les Fig. 8A a 8E montrent les deformations graduelles operees par un 
changement de coefficient d* anisotropic, sur un modele a trois facies genere 
par simulation sequentielle d' indicatrices; et 
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- les Fig. 9 A k 9E montrent les deformations graduelles locales d'un modele a 
trois facids, generees par simulation sequentielle d'indicatrices. 

Description detaiUee de la methode 

On considere une zone d' etude que Ton subdivise par une grille a N 
5 mailles et Ton cherche construire des realisations ou representations d'un 
modele stochastique d'une certaine grandeur physique z representant par 
example la permeabilite des formations dans la zone. Le modele recherche doit 
s'ajuster a des donnees obtenues par des mesures ou observations en un certain 
nombre de points, et notamment a des donnees non lineaires. 

10 Le calage d'un modele stochastique a des donnees non lineaires peut 

etre formulee comme un probleme d' optimisation. On designe par f^^= (fi*"^^ , 
p^abs ^ f^obs fp^*^*)^ le vecteur des donnees non lineaires observees ou mesur^es 

dans le milieu etudi^ (la zone reservoir), et par f = (fi , f^, fp), le 

vecteur correspondant des reponses du modele stochastique de la permeability 

15 Z = (Zi, Z2,..., Zn). Le probleme de contraindre le modele stochastique Z par 
les observations consiste h g^nerer une realisation z de Z qui reduit a une 
valeur assez faible une fonction objectif que Ton d^flnit comme la somme des 
erreurs quadratiques ponderees des reponses du modele par rapport aux 
observations ou mesures dans la zone-reservoir, soit : 

^ /-I 

ou coi reprdsente le poids attribue a la reponse f^ Les fonctions fi (i=l, 2, 
p) ainsi que la fonction objectif O sont des fonctions du vecteur Z. On se 
U"Ouve devant un probleme d*optimisation de dimension N. 
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Pour etendre le formalisirie developp^ dans le brevet deja c\t6, i la 
deformation graduelle de realisations gdndrees par simulation sequentielle non 
ndcessairement gaussienne, on va partir d'un vecteur gaussien 21 N variables Yi 
avec i = 1, 2, N, mutuellement independantes, de moyenne nulle et de 
5 variance egale ^ T unite, et on va definir N variables uniformes mutuellement 
independantes Ui, U^, U3,...Un par : ; 

Ui = G(Yi) V i=l,2, ...N 

oil G designe la fonction de repartition gaussienne standardis^e. 

Ceci pose, on va appliquer Talgorithme de deformation graduelle 

10 developpe dans un cadre gaussien, au vecteiu- gaussien Y = (Yi, Y2, , Yn), 

pour construire une chame continue de realisations du vecteur uniforme U = 

(Ui, U2, Un). Etant donne deux realisations independantes yaCt yb de Y, on 

definit la chafne de realisations u(t) du vecteur U que Ton obtient selon la 
relation : 

15 «(0 = G(y^cosr + .y^sinO (2) 

Pour chaque t, u(t) est une realisation du vecteur U. Ensuite. par 

echantillonnage de la distribution conditionnelle /;(zJCi=l. 2 N) en 

utiiisant les elements du vecteur u(t), on obtient un vecteur z(t) qui est, pour 
chaque t, une realisation du vecteur aieatoire Z. En consequence, on peut 
20 ajuster le parametre t conune dans le cas gaussien de fa?on caler 2(t) h des 
donnees non lineaires, Cette procedure est iteree jusqu'a atteindre un calage 
satisfaisant. 
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Calage d'on modele de facies a des donnees de pression obtenues par 
des tests de puits 

Pour illustrer I' application de la methode d' optimisation stochastique 
ci-dessus definie, on va chercher a caler un modele stochastique de reservoir li 
5 des donnees de pressions obtenues par des tests de puits. La construction du 
modele de reservoir est derivee d'une formation petroliere reelle qui comporte 
trois types de facies : deux facids reservoir de bonne qualite (facias 1 et 2) et 
un type de facies reservoir de trds mauvaise qualite (facias 3). La table 1 
defmit les proprietes petrophysiques des trois facies : 















facids 1 


10. 


10. 


10. 


17 


2.1857 


facias 2 


1. 


1. 


I. 


14 


2.0003 


facias 3 


0.1 


0.1 


0.001 


9 


1.8148 



Pour reprdsenter la distribution specifique de facies de la formation 
petrolifere, on genere d'abord une realisation binaire pour representer le facies 
3 et son complement. Ensuite, dans la partie complementaire du facias 3, on 
gendre une autre realisation binaire independante de la premiere pour 
15 representer les facias I et 2. La formation est discretis^e par un maillage 
regulier de 60x59x15 blocs de taille 15mxl5mxl.5m. Un modele de 
variogramme exponentiel est utilise pour estimer les distributions 
conditionnelles. La direction d'anisotropie principale est diagonale par rapport 
au maillage. Les portees du variogramme du facies 3 selon les trois directions 
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d'anisotropie sont respectivement de 300m, 80m et 3ra. Les portees du 
variogramme des facies 1 et 2 selon les trois directions d' anisotropic sont 
respectivement de 150m, 40m et 1.5m. Les proportions des facias I, 2, 3 sont 
respectivement de 6%, 16% et 78%. 

5 On a realise un test de puits au moyen d'un slmulateur de test de putts i 

differences finies, de la mani^re decrite par : 

Blanc, G. et al : « Building Geostatistical Models Constrained by Dynamic 
Data • A Posteriori Constraints » in SPE 35478, Proc. NPF/SPE European 3D 
Reservoir Modelling Conference, Stavenger, Norway, 1996, 

10 On voit sur la Fig.l la couche mediane d'une realisation qui sert de 

module de reference pour cette validation. La section du puits qui a 6i€ 
perforee traverse horizontaleraent la couche mediane du modele de reserv^oir le 
long de Taxe x. Le diamdtre du puits est de 7.85cm, la capacite du puits est 
nulle et les facteurs parietaux (skin) des facies 1 , 2 et 3 sont respectivement 0, 

15 3 et 50. Le test de puits synthetique dure 240 jours avec un debit constant de 5 
m3/jour de fagon a invesdguer presque la totalite du champ petrolifere. La 
Fig. 2 montre la variation de pression en fonction du temps. 

Le but ^tait de construire des realisations du modele de facias contraint 
par les facias rencontre le long du puits et par la courbe de pression obtenue 

20 lors du test de puits. La fonction objectif est d^fmie conmne la somme des 
differences quadratiques entre les reponses en pression du modele de reference 
et les reponses en pression de la realisation. Du fait que le componement 
dynamique du modele de reservoir est principalement contrdle par le contraste 
entre les facies reservoir de bonne et de mauvaise qualitds, la realisation 

25 binaire utilisee pour gdn^rer les facies I et 2 a ete fixee d'abord et seule la 
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realisation binaire utilisee pour genefer le faciSs 3 a €t6 deform^e pour le 
calage des donnees de pression. 

Les reponses en pression rdsultant des tests de puits pour les cinq 
realisations des Fig.3A i 3E, sont diferentes de celle du modele de reference, 

5 conime le montrent les Fig. 4 A £ 4E. Partant respectivement de ces 5 
realisations independantes, en utilis^t la methpde de calage iteratif ci-dessus, 
on obtient apr^s quelques iterations, cinq realisations calees (Fig. 5A a 5E) 
pour lesquelles les courbes de pression correspondantes sont tout a fait en 
accord avec celles du modele de reference, comme le montrent les Fig.6A h 

10 6E. 

Deformation graduelle par rapport aux parametres structuraux 

Dans de nombreux cas, on ne dispose pas de donndes suffisantes pour 
en deduire les parametres structuraux du modele stochastique : moyenne, 
variance, fonction de covariance, etc. Ces paramdtres structuraux sont souvent 

15 donnes en termes d'intervalles ou de distributions a priori. Si leurs valeurs sont 
fausses, il est vain de rechercher une^ realisation calee sur des donnees non 
lineaires. Aussi, il est essentiel pour les applications de pouvoir proceder a une 
deformation graduelle d'une realisation avec modification simultanee de 
nombres aleatoires et des parametres structuraux. L'ulgorithme de simulation 

20 sequentielle defmie par 1' Equation (I) foumit cette possibility de changer 
simultanement ou separdment Toperateur structurel S et le vecteur uniforme U. 
Les Fig.8A-8E montrent les deformations graduelles que Ton obtient pour une 
realisation fixe du vecteur uniforme U quand on change le coefficient 
d' anisotropie. 
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Deformation graduelie locale ou regionalisee 

Quand les observations sont repaities sur des zones differentes d*une 
formation etudiee, un calage qui utiliserait une deformation globale risquerait 
d'etre inefficace car I' amelioration de T accord que Ton obtiendrait sur une 
zone pourrait le det^riorer sur une autre zone. II est done bien preferable 
d'appliquer des deformations graduelles zone par zone. On consid^re une 
subdivision du vecteur U en un certain nombre n de sous-vecteurs \j\ 



U" mutuellement independants ce qui rend possible de realiser 
individuellement leur deformation graduelie. En appliquant separement i 

10 chacun des sous-vecteurs U\ U* U^ Talgorithme de deformation 

graduelie, on obtient une fonction de dimension n, du vecteur uniforme U : 



G(y>osf, + rJsinfj' 
GcK„^cosf^ + r/^sin t^) 



G(K;cosr„ + i;"sinfJ 



oil et Yl pour tout i =1, 2, n, sont des sous-vecteurs gaussiens 
independants. Pour un ensemble de realisations de et F/, on resout un 
15 probleme d' optimisation de n paramdtres ti, t^, tn pour obienir une 
realisation qui ameliore, sinon maintient, le calage des donnees. Cette 
procedure peut etre itdree jusqu'i atteindre un calage satisfaisant. 

En realisant des deformations locales graduelles, on ameliore de fa(;on 
significative la vitesse de calage dans tous les cas ou les mesures ou 
20 obser\'ations sont reparties dans differentes zones du milieu. 
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On voit bien Teffet de cette deformation locale graduelle sur le module 
i trois facias, des Fig. 9A i 9E ou seule la partie inferieure gauche delimitee 
est affectee. 

La methode selon T invention pent etre facilement generalisde h la 
5 deformation graduelle d'une representation ou realisation d un modele 
stochastique quelconque, car gdnerer une realisation d'un tel modele 
stochastique se ram^ne toujours a la generation des nombres uni formes. 
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REVENDICATIONS 



1) Methode pour deformer graduellement une representation ou 
5 realisation, gdner^e par simulation sequentielle, d'un modele stochastique non 

necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre i un ensemble de donnees 
collectees dans le milieu par des raesures et observations prealables, relatives h 
son etat ou sa structure, caracterisde en ce qu'elle comporte T application d'un 

10 algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) ^ N variables mutuellement independantes qui est relie 5 un 
vecteur uniforme U h N variables uniformes mutuellement independantes par 
la fonction de repartition gaussienne (G), de fagon h construire une chaine de 
realisations u(t) du vecteur U, et T utilisation de ces realisations u(t) pour 

15 generer des realisations z(t) de cette grandeur physique, que Ton cale par 
rapport aux donnees. 

2) Methode selon la revendication 1, caracterisde en ce que Ton definit 
une chaine de realisations u(t) du vecteur (U) ^ partir d'une combinaison 
lineaire de realisations du vecteur gaussien (Y) dont les coefficients de 

20 combinaison sont tels que la somme de leurs carres est egale a T unite. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle du modele representatif du milieu 
heterogdne, simultanement par rapport aux parametres structuraux et aux 
nombres aleatoires. 
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4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle separee d'un nombre n de parties du 
modele representatif du milieu heterog^ne tout en preservant la continuite 
entre ces n parties du meddle, par une subdivision du vecteur uniforme en n 
sous-vecteurs mutuellement independants. 
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Fig.5A Fig.5B Fig.SC Fig.5D Fig.SE 
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